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ABSTRACT 

The high-rcdshift Universe contains luminous Lya emitting sources such as galaxies 
and quasars. The emitted Lya radiation is often scattered by surrounding neutral 
hydrogen atoms. We show that the scattered Lya radiation obtains a high level of po- 
larization for a wide range of likely environments of high-redshift galaxies. For example, 
the back-scattered Lya flux observed from galaxies surrounded by a superwind-driven 
outflow may reach a fractional polarization as high as ~ 40%. Equal levels of polar- 
ization may be observed from neutral collapsing protogalaxies. Resonant scattering 
in the diffuse intergalactic medium typically results in a lower polarization amplitude 
( <7%), which depends on the flux of the ionizing background. Spectral polarimetry 
can differentiate between Lya scattering off infalling gas and outflowing gas; for an 
outflow the polarization should increase towards longer wavelengths while for infall the 
opposite is true. Our numerical results suggest that Lya polarimetry is feasible with 
existing instruments, and may provide a new diagnostic of the distribution and kine- 
matics of neutral hydrogen around high-redshift galaxies. Moreover, polarimetry may 
help suppress infrared lines originating in the Earth's atmosphere, and thus improve 
the sensitivity of ground-based observations to high-redshift Lya emitting galaxies 
outside the currently available redshift windows. 
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exploit this window further by pushi ng to even higher 
redshifts and fain ter flux levels (e.g. IStark et al.l l2007bl : 
iNilsson et al.ll2007l . and references therein). 

The interpretation of existing and future observations 
is complicated by the fact that Lya photons are typically 
scattered both in the interstellar medium (ISM, e.g. Hansen 
fc Oh 2006) as well as in the intergalactic medium (IGM, 
Loeb fc Rybicki 1999; Santos 2004; Dijkstra et al 2007b). 
The impact of scattering on the observed Lya flux and 
spectrum may be derived by careful modeling of the ob- 
served Lya line proflle, combined with constraints on the 
gala xy's stellar population derive d from its broad band col- 
ors l|Verhamme et al"]|2006l . [20081 ). Also, observations of the 
Lya line in local starburst galaxies can shed light on the 
processes that regulat e the transfer and escape of Lya from 
star forming galaxies (|Kunth et al.lll998l ; lHaves et al.ll2007l ). 
Of course, it is not obvious whether local starburst galax- 
ies are representative of high-redshift Lya emitters (LAEs), 
and detailed modeling of the line profile may be hampered 
by the quality of the data of high-redshift LAEs. 

In this paper, we describe a complementary approach 
to studying the environments of high-redshift LAEs, which 
involves the polarization properties of the scattered Lya 
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radiation. In the context of the solar system, the po- 
larization of resonantly scattered Lya photons from the 
sun by neutral hydrogen in the interplanetary medium 
and in the Earth's atmosphere has bee n studied exten- 
sively (e.g. iBrandt fc ChamberlainI Il959l ). occasionally by 
using a Monte -Carlo method to c ompute the Lyg r adia- 
tive transfer (iModali et all 1 19721: [K eller et ajj 119811 ). In 
the cosmological context. iLee fc Ahn. ( 1998 ) showed that 
Lya emerging from an unresolved star-bursting galaxy 
may be polariz e d to a considerable level (~ 5%), and 
iRvbicki fc Loebl l|l999l) showed that scattering of Lya 
photons from a point source embedded in the Hubble 
flow of a neutral IGM would produce a polarized Lya 
halo around the source. Despite the emergence of several 
Monte-Carlo Lya radiative transfer cod es in recent years 
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polarization properties of scattered Lya photons around 
galaxies have not been investigated in a broader context 
beyond this early work. 

The goal of this paper is to demonstrate that Lya ra- 
diation may be highly polarized around high-redshift galax- 
ies in a broad range of cosmological circumstances. We will 
show that Lya polarimetry may place constraints on the 
kinematics of the gas through which the photons propagate. 
This is important: understanding the impact of scattering 
in both the ISM and IGM on the observed Lya properties is 
required to fully exploit Lya observations as a probe of the 
high-redshift Universe. Furthermore, we discuss the possibil- 
ity of separating Lya sources from low-redshift line emitters 
(such as [on], [OIII], Ha, H/3 emitters, and OH-skylines) 
based on polarimetry. 

The outline of this paper is as follows: In §[5] we describe 
the basic principles of Lya radiative transfer and polariza- 
tion, and how polarization is incorporated in our radiative 
transfer calculations. In § [3] we present our main numeri- 
cal results. In § |4] we discuss our results and their impli- 
cations, before presenting our final conclusions in § [5] The 
parameters for the background cosmology us ed throughout 
this p aper are {Q^Ma, h) = (0.24, 0.76, 0.73) l|Spergel et all 
l2007l ). 



2 LYa SCATTERING & POLARIZATION 

2.1 Lya Radiative Transfer Basics 

We start by summarizing the basic principles of Lya scatter- 
ing. We express photon frequency u in terms of a dimension- 
less variable x = {v ~ vq)/ Avd, where Avd = VthVo/c, and 
Vth is the thermal velocity of the hydrogen atoms in the gas, 
given by vth = sJ^ksT jm^, where fcs is the Boltzmann con- 
stant, T — 10* K the gas temperature, nip the proton mass 
and vq = 2.47 x 10^^ Hz is the Lya resonance frequency. For 
reference, the optical depth through a column of hydrogen, 
A'hi, for a photon in the line center, ro, is given by 

iVm \{ T \ -0-5 



The optical depth for a photon at a frequency x reduces to 
dy 

T"0 ^ J-oo (j/ - + °^ 



/C 
-c 



(e.g. iRvbicki fc LightmanI 1 19791 ) where a is the Voigt 
parameter given by a = A^ij ^"k^vb = 4.7 x 
10"*(«th/13 km s-^)-\ where A21 = 6.25 x 10* s"^ is the 
Einstein A-coefRcient for the transition. The transition be- 
tween 'wing' and 'core' occurs &% x = Xt k, 3.3 for T — 10* 
K. 

Under most astrophysical conditions, absorption of a 
Lya photon is followed be re-emission of a photon of the 
same energy in the frame of the atom. However, in the ob- 
server's frame the Lya photon's energy before and after scat- 
tering is modified by Doppler shift due to the thermal mo- 
tion of the atom, and scattering is only 'partially' coherent 
in the observer's frame. These Doppler shifts are important 
as they cause the photon's frequ ency to change by an rms 
shift of Ai/_D per scattering event l|Osterbrockll 1963 ) . There- 
fore, Lya scattering through an optically thick medium can 
be described as a sequence of rando m walks in both fre- 
quency and real space (e.g rHarringtor^ [l973l : lNeufeld|[l990l : 
iLoeb fc Rvbickilll999l ). Frequency diffusion is very efficient 
and spatial diffusion occurs predo minantly when photons are 
in th e wing of the line profile (e.g. lAdams|[l972l : iHarringtonl 
Il973h . Hence, the last scattering event occurs in the wing 
of the line profile for the majority of Lya photons emerging 
from an optically thick medium q 

Other quantities of relevance to our discussion are: it) 
the scattering phase function (also known as the anisotropy 
function) which gives the probability p(0) that the scattered 
photon is re-emitted at an angle d relative to the incoming 
photon; and (ii) The degree of polarization n(S) caused by 

^pji, where 7|| and 



scattering which is defined as n(0) e 

7j_ are the intensities parallel and perpendicular to the plane 
of scattering (defined by the wave vectors of the ingoing and 
outgoing photons), as a function scattering angle Q. Both 
p(6') and n(S) are different for core and wing scattering, as 
we discuss in more detail below. 



2.2 Phase Function &: Degree of Polarization for 
Lya Scattering 

For resonant scattering the phase function and the degree 
of polarization are strongly dependent on the atomic lev- 
els involved in the scattering event, and must be calculated 
quantum-mechanically. For example, the sequence of scat- 
tering events 15i/2 'iP\ii 15i/2 result in an un- 
polarized isotropically re-emitted Lya photon, while the 



scattering events 15*1/2 



2P. 



3/2 



15*1 /2 produc e Lyg 



core; 
wing, 



(2) 



with a maximum polarization of | (e.g. iHamiltonl 1194*71 : 
IChandrasekIiailll960l : I Ahn et al.|[200^ ). Here we used a no- 
tation nLj, in which n, L and J are the principal quantum 
number, orbital angular momentum number, and total an- 
gular momentum number of the hydrogen atom involved 
in the scattering event, respectively. By summing over all 
possible Lya transitions, the phase function and degree of 



^ Note that Lya photons typically scatter only ~ ro times before 
emerging from a medium of optical depth tq, as opposed to the 
TQ-scaling, that is expected in the absence of frequency diffusion. 
For scattering events that occur in the wing, a so called 'restoring 
force' pushes the photon b ack towards the core by an average 
amount of -1/a; llOsterbroc k 1962; Adams 1972). 
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Figure 1. Schematic geometry of the last-scattering event of a 
Lya photon that occurs at a radius vector r away from the source 
galaxy. Before last scattering the photon's propagation direction 
is n (which does not necessarily lie in the plane of the sky). After 
scattering, the photon travels in a direction n' which is perpen- 
dicular to the sky plane. The polarization vector e' after scatter- 
ing lies in the plane of the sky and intersects the projected (on 
the sky plane) radius vector r at an angle x- This photon con- 
tributes cos^ X siii^ X to the intensity of the radiation field 
parallel (/;(a)) and perpendicular {Ir{a)) to r, respectively. Here 
a = r/dAiz) (see text for additional details). 



polarization for resonant scatterinR were found to be (e.g. 
iBrandt fc Chamberlainlll959l : iBrasken fc Kvrolalll998l ) 



p(e) = \ COS 

' 12 12 



6, n{6) = 



f + cos2 e ' 



(3) 



As is shown in Appendix [A) this corresponds to a case of a 
superposition of Rayleigh scattering plus isotropic scatter- 
ing with corresponding weights of 1/3 and 2/3. The phase 
function p{9) satisfies J p{9)dfl — 4n. 

In quantum mechanics a Lya scattering event cannot 
go separately through either 15*1/2 — > 2Pi/2 — > 151/2 or 
15i/2 2P3/2 — > 15i/2. Instead, each scattering event is 
a sup erposition of sc attering through both levels simultane- 
ously. IStenfld (|1980D has shown that this introduces quan- 
tum interference terms which cause Lya wing-scattering to 
be described by Rayleigh scattering for which 



p{o) = l + leasee, u{e) 



(4) 



1 + cos2 e ■ 

This remarkable result implies that Lya wing scattering can 
produce th ree times mo re polarization than Lya resonant 
scattering (|Stenfldll980l 'l. 



2.3 Incorporating Lya Polarization in Radiative 
Transfer Codes 

The code we used for the Monte-Carlo L ya radiative trans- 
fer is described in IPiikstra et al.l l|2006t ). The code follows 
individual Lya photons through spherical concentric shells 
with user-specified density and velocity fields. The code ac- 
curately describes the process of frequency and spatial diffu- 
sion in Lya radiative transfer as described in § 12.11 (includ- 
ing other effects that were not discussed such as for exam- 
ple energy loss due to recoil). For a detailed description of 
the co de, the interested reader is referred to iDiikstra et al] 
The code was adapted to the context of this paper 
as follows: 
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Figure 2. Schematic illustration of the polarization of scattered 
Lycf radiation. Arrows pointing away (towards) the galaxy rep- 
resent a radiation field for which /j < (i.e. V{a) > 0), while 
II > Ir is represented by arrows tangential to spheres of constant 
radius. The magnitude of the polarization is represented by the 
size of the arrows. The arrows vanish for unpolarized radiation. 
In all cases considered in this paper, 7^(0) < (this corresponds 
to a Stokes parameter Q < 0), and the polarization vectors form 
concentric shells surrounding the central source. 



• The modified code uses different phase functions for 
scattering in the core and wing as described in § 12.21 The 
original code (which was constructed for a different purpose) 
assumed only the Rayleigh scattering phase function. 

• In § 12.21 we showed that the polarization of scattered 
Lya photons depends strongly on whether Lya was reso- 
nantly scattered or not. We also mentioned in § 12.11 that 
the transition between resonant and wing scattering occurs 
at X 3.3. However, it is possible to determine whether 
a photon scattered resonantly for each scattering event di- 
rectly from the Monte-Carlo simulation. While generating 
the velocity of the atom that scatters the Lya photon, a 
scattering event is defined to be resonant if it occurred less 
than Xcvit = 0.2 Doppler widths away from resonance m 
the frame of the atom involved in the scattering event (see 
Appendix |X] for a justification of this number). The precise 
choice of a;crit does not affect our main results. 

• The code was modified t o include polarization based 
on the scheme developed by iRvbicki fc Loebl (|l999'). Al- 
though their method is applicable strictly to wing/Rayleigh 
scattering, only minor modifications are required to calcu- 
late polarization for resonantly scattered Lya photons. This 
is mainly because resonant scattering is a superposition of 
isotropic and R ayleigh scattering. We will briefly discuss the 
iRvbicki fc Loeb (,1999. ) method, before describing our mod- 
ifications; 

IRvbicki fc Loebl (|l999l ) assigned 100% linear polarization 
to individual photons by endowing each Lya photon with a 
polarization vector e, which is perpendicular to the photon's 
propagation direction n, i.e. e • n = 0. In this formulation, 
the Stokes parameters result from binning together multiple 
independent photons. 

The propagation direction (n') of a photon after scatter- 
ing is obtained by means of a rejection method: A random 
direction n', and a random number R € [0, 1] are gen- 
erated. The new propagation direction is accepted when 
R < 1 — (e ■ n')'^ (as shown in Appendix IA3I this corre- 
sponds to a more general case of the phase-function p{0) 
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given in Eq|4ll. The polarization vector after scattering (e') 
is obtained by normalizing the vector g', which is obtained 
by projecting the polarization vector prior to scattering (e) 
onto the plane normal to n'. Symbolically, e' = g/|g|, where 
g = e — (e ■ n')n'. 

The observed radiation is characterized by its intensities 
parallel (/;) and perpendicular (Ir) to the radius vector 
to the location of last scattering, denoted by r, projected 
on the plane of the sky (see Figure [!}, denoted by f. For 
the last scattering event n' is perpendicular to the plane 
of the sky, and the polarization vector e' lies in the plane 
of the sky (see Fig [1} . We use x to denote the angle be- 
tween the polarization vector e' and f . A photon that is ob- 
served from impact parameter a = r/dA{z) (with dA(z) be- 
ing the angular diameter distance to redshift z) contributes 
cos^ X to /; and sin^ x t o Ir for Rayleigh (or wing) scattering 
l|Rvbicki fc Loebl ll999V 

However, for resonant scattering a photon contributes 
Prcos^X + |(1 — Ph) to /; and pRsin^x + |(1 — Pr) to 
Ir- Here, pa = {l + cos^9)/{^ + cos^ 9) is the probability 
that a core photon was Rayleigh scattered in the last scat- 
tering event (see Appendix Al), in which 6 is the angle be- 
tween the propagation directions of the incoming and outgo- 
ing photons at last scattering (cos 6 = n ■ n'). The fact that 
PR depends on angle 6 has a simple reason. Although reso- 
nant scattering may be viewed as a superposition of Rayleigh 
scattering and isotropic scattering of weights 1/3 and 2/3, 
these weights are averaged over solid angle. The actual 
weight of each process is angle dependent (see Appendix IX)) . 
In other words, in the case of resonant scattering the intensi- 
ties in radial and azimuthal directions are obtained by scal- 
ing down the intensities obtained from Rayleigh scattering 
by a factor of pR, and by adding an unpolarized intensity of 
magnitude 1 — pR. 

The observed fractional polarization at impact parameter 
a is given by 



P(a) = 



lilies) 



Ii{a)+Ir{a) 



(5) 



The fractional polarization can also be expressed in terms 
of the Stokes parameters Q and I as V — \Q\/I- In this 
paper, we always find that Ji < Ir, and Q < 0. This may be 
represented by arrows drawn in concentric rings around the 
central Lya emitter, where the size of an arrow indicates the 
magnitude of V (see Fig[2ll. 

• Monte-Carlo calculations of Lya transfer can be accel- 
erated by skipping scattering events that occur in the line 
core (e.g. lAhn et al.i[2002l : iDiikstra et al.ll200d '). One has to 
be somewhat careful when applying this technique in calcu- 
lations of polarization, as it reduces the number of scatter- 
ing events that a Lya photon encounters and hence distorts 
its resulting polarization state. However, we have verified 
explicitly that the polarization properties of Lya emerg- 
ing from optically thick clouds are correctly calculated by 
the accelerated scheme. Generally speaking, the accelerated 
Monte-Carlo scheme yields the correct polarization when- 
ever this scheme reliably computes other observable quanti- 
ties such as the spectrum and surface brightness profile. 

Lastly, we point out that the expressions for p{9) and 
Tl{9) in § 12.21 were derived under the assumption that the 
radiation field prior to the scattering event was unpolarized. 



The precise functional forms of p{9) and H(6') are differ- 
ent when the radiation field prior to scattering is polarized 
(see Appendix Al). The procedure of propagating the po- 
larization vector in a single Rayleigh scattering event that 
was outlined above, naturally accounts for the polarization 
dependence of Ii{9) and p{9) (see Appendix A3). For res- 
onantly scattered Lya the situation is more co mplex and 
the density matr i x formalism shou ld be used (e.g. lLeelll994l : 
iLee fc Ahnlll998l : lAhn et al.ll2002l ). However, the precise de- 
pendence of p{9) and H(S) on initial polarization does not 
affect the calculations regarding resonant scattering. The 
maximum fractional polarization we find for resonantly scat- 
tered Lya is 7% (see §[3]2}. The typical polarization of the 
Lya radiation field 'seen' by atoms involved in the scattering 
process is even lower. This only introduces changes in the 
phase function at the level of a few per cent, and we have 
verified that our results regarding resonantly scattered Lya 
are insensitive to the precise form of the phase function. 



3 LYa POLARIZATION AROUND 
HIGH-REDSHIFT SOURCES 

Next we calculate the expected polarization of the Lya line 
in a set of models which span the likely environments of 
high-redshift galaxies. Unless otherwise stated, we assume 
the redshift of Lya sources to be z = 5.7 in all cases. The 
Monte-Carlo radiative transfer calculations rely on stacking 
individual photons into bins, and the errorbars shown in 
some figures were calculated assuming Poisson fluctuations 
in the number of photons within a given bin. 



3.1 Backscattering off Galactic Superwinds 

3.1.1 Description of the Model 

A major fraction of Lyman break galaxies (LBGs) at high 
redshifts sho w evidence of b eing surrounded by outflowing 
enriched gas iShaplev et al.| [2003'). Scattering of Lya pho- 
tons by neutral hydrogen atoms in these outflows may cause 
the observed Lya line to be redshifted systematically r elative 
to the systemic velocity of the galaxy l|Ahn et al.ll2003l ) . This 
redshift is attributed to the Doppler boost that Lya photons 
experience when they scatter off the outffow on the far side 
of the galaxy back towards the observer (hence the term 
'backscattering'). This redshifted Lya flux is less prone to 
resonant scattering in the IGM, and is therefore more easily 
observable than Lya photons that were not backscattered. 

The impact of the outflow on the observable proper- 
ties of the Lya line depends on various parameters includ- 
ing the outflow speed, Wexp, the total column density of 
neutral hydrogen atoms in the outflow, A^hi, and the dust 
content and distribution insi de the outflow (|Hansen fc "Ohl 
I2OO6I : IVerhamme et al.l I2OO6I ) . We show results for a sin- 
gle expanding thin shell of gas with A^hi ~ 10^^ cm~'^ 
and A'hi = 10'^° cm~^, based on col umn densities that 
are observed in Lya e mitting galaxies (|Kunth et al.l [l998l : 
IVerhamme et al1l2008h . We discuss the choice of our model 
in more detail after presenting our results below. The out- 
flow speed is typically of order t he circular velocity of th e 
host dark matter halo, Vdrc (e.g. iFurlanetto fc Loebll2003l ). 
We assume the galaxy to be embedded in a halo of total 
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8160 
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Frequency (x) 



Figure 3. The observable properties of the Lya radiation emitted by a galaxy which is surrounded by a single thin spherical shell of 
gas that expands at a speed Voxp = 200 km and contains neutral hydrogen with a column density of either A^jji = 10^^ cm"'^ (red 
dashed lines) or Afjji = 10'^'' cm~^ {solid lines). These parameter choices represent the generic conditions in galaxies surrounded by a 
superwind driven outflow (see text). The dotted line in the left panel shows the Gaussian Lyo emission line prior to scattering. The figure 
shows that backscattering off the expanding outflow results in a systematic redshift of the Lyo line. The top and bottom right panels 
show the surface brightness proflle (units are erg s"'^ cm~^ arcsec"^) and fractional polarization 'P{ct) as a function of a, respectively. 
Here, the amplitude of a and the surface brightness S{a) are written in terms of etc, which is related to the distance of the expanding 
shell from the galaxy (see text). The Figure shows that the polarization reaches values as high as ^ 40%. 



Njjj=10^^ cm ^ 




5 10 
Number of scattering events 



15 



Figure 4. The histogram shows the angle-averaged polariza- 
tion as a function of the number of times Ly« photons scatter 



for the superwind model with A^hi 



10^ 



cm ^. The angle- 



averaged polarization is largest for photons that scatter only once 
((P(a)) ~ 30%), and almost vanishes for A^scat ^ 10. The thick 
solid line shows the probability distribution of the number of 
times Lya photons scatter. Clearly, the majority of photons scat- 
ter only once, which causes the observed Lyo flux to be highly 
polarized. 



mass Mtot 
km s~^. 



3 X 10 Mq, which corresponds to Vdi 



200 



X (A) 
8140 8160 



8180 




-50 
Frequency (x) 



100 



Figure 5. The frequency dependence of the polarization at « = 
{thin black histogram), 0.5ac {red dashed histogram), and etc {blue 
dotted histogram). The polarization increases with impact param- 
eter. The thick solid line shows the observed spectrum which was 
also shown in Fig |3] (the units are arbitrary). The bulk of the 
photons have frequencies —50 < x < 10. Within this frequency 
range, the polarization increases towards the red. This is because 
the redder Lyo photons appear farther from resonance in the 
frame of the expanding gas; consequently they scatter less and 
achieve a higher polarization amplitude (see Fig|4]l. This generic 
frequency dependence of the polarization amplitude can be used 
as a fingerprint of outflows. The frequency dependence is reversed 
for infalling gas (see § I3.3|l . 



3.1.2 Results 

The emitted Lya spectrum prior to scatterin g is assumed to 
be a G aussian with a velocity width a = Vdrc l|Diikstra et al.l 
l2007bl ) and is denoted by the dotted line in the left panel 
of Figure |3l The dashed (solid) lines represent the observed 
spectrum for an outflow with Nni = 10^^ cm"'^ (10'^° cm~^). 
The Figure shows that backscattering off the expanding out- 



flow results in a systematic redshift of the Lya line (as 
mentioned above), with an increase in redshift as A^hi in- 
creases. The top right panel shows that the surface bright- 
ness profile remains flat out to some value oi a — ac af- 
ter which it drops to zero (as there is no scattering out- 
side the shell). The magnitude of ac and the surface bright- 
ness S{a) = Ii{a) + Ii{a) are determined by the radius 
of the expanding shell. For example, if the outflow is lo- 
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cated at a distance d = 10 kpc from the central galaxy 
(and its thickness ^ d), then the surface brightness pro- 
file drops to zero at Qc — 1.6"(d/10 kpc) for a galaxy at 
z — 5.7. The mean surface brightness scales approximateljQ 
as {S)^Xx 10-1** (Qc/1.6")-^(iLy«/10''^ erg s'^) erg s"! 
cm~^ arcsec"'^, where X = 3 for A''hi ~ 10^^ cm~^, and 
A" = 5 for iVm = 10^" cm'^ 

The bottom right panel shows that the fractional po- 
larization, T'{a), increases towards larger values of ct and 
reaches a maximum near a^. The polarization amplitude 



10^ 



and 



obtains values as high as ~ 40% for A*'hi 
~ 18% for Nm = 10^" cm'^ 

The reason for the overall large values of the polariza- 
tion is easy to understand. Although the line center optical 
depth through the outflow is very large for both models, 
TO = 5.9 X lQ^{Nm/lQ^'^ cm"^) (see Eq[T|, most photons 
reside in the wings of the line profile when they reach the 
expanding shell. For example, in the model with A'hi = 10^^ 
cm~^, the optical depth is less than unity for 63% of the pho- 
tons when they enter the outflow for the flrst time. There- 
fore, most photons scatter only once in the wing of the line 
profile before escaping to the observer. 

The above interpretation is quantified by Figure |4l in 
which the thick solid line shows the probability distribution 
for the number of scattering events encountered by the Lya 
photons, P(iVscat). The function P(iVgcat) peaks at A'^^cat = 1 
and rapidly decreases for increasing A'scat- Also shown as 
the histogram is the angle- averaged polarization of the pho- 
tons as a function of TVscat. The polarization is largest for 
photons that scatter only once {{P{a)) ~ 30%), and nearly 
vanishes for A'scat ^ 10. The decrease in the polarization am- 
phtude with increasing Ascat reflects that the radiation fleld 
becomes increasingly isotropic as the number of scattering 
events increases. As Ahi increases, the fraction of photons 
that escape after a single scatter decreases, and the overall 
polarization declines. 

Figure \5\ shows the frequency dependence of the polar- 
ization at three different impact parameters, namely a — 
{thin black histogram), 0.5ac {red dashed histogram), and Qc 
{blue dotted histogram). As already illustrated in Figure |3l 
the polarization increases with increasing impact parame- 
ter. The thick solid line shows the observed spectrum which 
was previously shown in Fig [3] The bulk of photons have 
frequencies —50 < a; < 10. Within this frequency range, the 
polarization increases towards the red. This is because the 
redder Lya photons appear farther from resonance in the 
frame of the gas; they therefore scatter less and achieve a 
higher polarization amplitude (see Fig|4]). This frequency 
dependence of the polarization can be used as a flngerprint 
to distinguish outflows from infall (see §[33J. 

One may wonder whether the above results depend sen- 
sitively on the assumed model: e.g. do the result change 
for a continuous, dusty, wind with a radial dependence of 
its outflow speed? A wind for which the outflow speed in- 
creases linearly with radius w ould resemble th e IGM with 
a Hubble flow as discussed in iRvbicki fc Loebl l|l999D . who 
found even higher levels of polarization. Furthermore, scat- 



^ Fig \3\ shows that the surface brightness remains almost con- 
stant out to a/ac ~ 0.8. The mean surface brightness, (5), was 
therefore estimated by dividing the total scattered flux by 7ra^. 



tering through an optically thick collapsing gas cloud also 
results in a comparable polarization (Pmax ~ 35%, see ij |3.3p . 
Changing the sign of the velocity of the gas only affects the 
frequency dependence of the polarization (see § 13. 3p . The 
column density of neutral hydrogen in both these models is 
A'hi ^ 10^" cm~^, which implies that polarized Lya is also 
expected for outflow models with larger column densities 
and non-zero velocity gradients. 

Furthermore, the presence of dust in the expanding 
wind can boost the polarization. The reason is that pho- 
tons that scatter only once obtain the highest level of polar- 
ization, while it decreases for photons that scatter multiple 
times (Fig|4]). However, photons that scatter multiple times 
traverse a longer path through the wind, which enhances the 
probability for absorption by dust. Hence, dust can prefer- 
entially quench the low polarization photons, which would 
strengthen our results. 



3.2 Resonant Scattering in the IGM 

3.2.1 Description of the Model 

Even after reionization, residual neutral hydrogen in the 
IGM can scatter up to ^90% of the Lvq emitted by galax- 
ies ou t of our line of sight jSantosI 20041 : iDiikstra et al.l 
l2007bl ). This scattered Lya radiation would appear as a 
faint Lyg halo surroun ding the central Lya emitting galaxy 
jhoeh fc Rvbickil 19991 ). We note that back-scattered Lya ra- 
diation is much less prone to resonant scattering in the IGM. 
Therefore we now consider galaxies that are not surrounded 
by wind-driven outflows, and examine examples that com- 
plement the cases discussed previously. 

We adopt the fiducial model described bv lDiikstra et al.l 
(|2007bl ) to examine the impact of the IGM on the observed 
Lya line profile of a 2; = 5.7 galaxy embedded in a dark 
matter halo with a total mass of Mtot = 10^^ Mq. The halo 
has a virial radius Tvir = 23 kpc and a circular velocity 
'^circ ~ 133 km s~i. The galaxy is assumed to form stars 
at a rate of M, — WMq yr~^ , and the escape fraction of 
its ionizing photons is taken to be /esc = 0.1. The result- 
ing LyQ luminosity of the galaxy is Li^ya— 2 x 10*^ erg s"^ 
for a gas metallicity oi Z = O.O5Z0. Gas within the virial 
radius is assumed to be fully ionized {x-n — 0) and no scat- 
tering occurs here. Just outside the virial radius the IGM 
density, p(rvir) ~ 20p, where p is the mean baryon density 
of the Universe. The density decreases as p(r) oc r~^ until it 
approaches p. Just outside the virial radius, the IGM is as- 
sumed to be collapsing at vigm ~ Vdic. Only at sufficiently 
large radii (r > lOrvir) the IGM is comoving with the Hub- 
ble flow (see Barkana 2004 and Dijkstra et al 2007b for a 
more detailed description and a motivation of this modeljf]. 

We also examine a case where we artificially boost the 



^ iDiikstra et al. I ll2007bl) account for gas clumping in the IGM. 
Most of the volume of the IGM consists of interclump gas that is 
underdense. For this reason the mean IGM opacity to Lyo pho- 
tons is lower in a clumpy IGM at a fixed value of the ionizing 
background (Dijkstra et al 2007b). Because clumps are not in- 
cluded in the radiative transfer code, we increased the ionizing 
background (compared to the clumpy case) in order to match the 
obse rved mean transmission of the IGM at a redshift z, (e~^)(z) 
(e.g. iFan et al.ll2006i) . 
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Figure 6. The observable properties of Lyo radiation that is resonantly scattered in the intergalactic medium (IGM). The cases 
considered here represent galaxies that are not surrounded by a superwind-driven outflow and complement the examples shown in Fig|3] 
(see text). The IGM is assumed to have a peculiar infall velocity (see text). The upper left panel shows the spectrum of Lya radiation 
that was transmitted without scattering once as the red dashed line. This spectrum reflects the line-emission spectrum observed from 
the galaxy (see text). The spectrum of resonantly scattered Lya radiation is shown as the solid line (both spectra are normalized). The 
upper right panels show that the fractional polarization is only 7-" (a) ~ 2%. The main reason is that each photon resonantly scatters 
~ 50 times, as shown by the dotted line. The lower panels show the results for the same calculation but with an increased level of the 
ionizing background. The boost in the ionizing background increases the transmission on the blue-side of the line {lower left panel). In 
this case each Lya photon scatters on average only a few times, and the scattered Lya is polarized up to 7% (lower right panels). 
Thus, the polarization of resonantly scattered Lyo photons in the IGM is low compared to that in § 13.11 and its level depends on the 
local level of the ionizing background. 



photoionization rate in the IGM. This boost may represent: 
(i) an enhanced local ionizing backgr ound due to clustering 
of undetected surrou nding sources (|Wvithe fc Loebl boosl : 
iDiikstra et aLll2007bl ): (ii) an enhanced ionizing background 
due to the presence of a nearby quasar; (in) a galaxy at a 
lower redshift where the residual neutral hydrogen fraction 
in the IGM is lower; (iv) an enhanced local photoionization 
rate due to vigorous star formation in the galaxy itself, or 
due to an enhanced escape fraction of ionizing photons; [v) 
a lower neutral fraction because of a lower overall density of 
hydrogen nuclei along the line of sight (e.g. for a galaxy on 
the edge of a void.) 

3.2.2 Results 

In the top left panel of Figure |6] we show the spectrum of 
Lya radiation that was transmitted without scattering once 
[dashed line), which represents the Lya spectrum observed 
from the galaxy. This spec trum was the f o cus of the analy- 
sis of e.g Santos (2004) and lDiikstra et all (|2007bl ). The plot 
shows that the blue side of the Lya line observed from the 
galaxy (x > 0) is eliminated by the IGM. The suppression 
extends into the red part of the line because of resonant scat- 
tering by residual neutral hydrogen gas that is falling onto 
the galaxy. Therefore, the peak of the observed spectrum 
is redshifted relative to the true line center by an amount 



which is set by the gas infall velocity. The redshift in veloc- 
ity is i;iGM(T'vir) = —Vciic. The spectrum of the resonantly 
scattered Lya halo is shown as the solid line. 

The top right panel of Figure|^shows the surface bright- 
ness profile of the scattered Lya halo as a function of a. The 
observed surface brightness of the scattered Lya photons is 
~ ]^g-i9 gj-g g-i cm~^ arcsec"^ near the center, and de- 
creases with increasing angular separation from the galaxy. 
Overall, the fractional polarization is small, ■P(a) ~ 2%. 
Besides the fact that resonantly scattered Lya radiation is 
always less polarized than the radiation in the line wings (see 
ij |2.2p . there is the fact that when the Lya enters resonance, 
it typically scatters multiple times before escaping towards 
the observer. This is shown by the dotted line which gives 
the mean number of scattering events for the Lya photons, 
{Nsca.t), as a function of impact parameter a. In this model 
(A'^scat) ~ 50. Hence, the Lya radiation field of resonantly 
scattered Lya photons will locally be close to isotropic, re- 
sulting in a low polarization. 

The lower panels in Figure |S] show the results for the 
model in which the ionizing background was boosted. The 
lower left panel shows that boosting the photoionization 
rate decreases the fraction of residual neutral hydrogen gas, 
which in turn increases the transmission on the blue-side 
of the line. The lower right panels show that the surface 
brightness profile falls slightly faster with a slight peak in 
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Figure 7. The observable properties of hya radiation emerging from a neutral collapsing gas cloud. The model represents protogalaxies 
in the process of their assembly. The left panel shows that the spectrum emerging from such a cloud has a systematic blueshift. The right 
panel shows that the Lya surface brightness profile (units are erg s~^ cm~^ arcsec"^) is rather fiat {top), and that it is highly polarized. 
The fractional polarization increases roughly linearly towards the edge of the cloud, and reaches a maximum value of Pmax ~ 35%. 



the surface brightness profile near the virial radius. The rea- 
son for this peak is that in our model the IGM is densest 
right outside the virial radius, which results in the largest 
effective optical depth here (this also produces the dip in 
the spectrum). In this model each Lya photon scatters on 
average only a few times which causes the scattered Lya to 
be polarized up to ~ 7%. The mean polarization decreases 
at Q < CYvir ~ 4". This is because in our model no scatter- 
ing occurs inside the virial radius, which extends out to a^ir- 
Therefore, the relative number of photons that was scattered 
at almost-right-angles decreases at a < Qvir, which in turn 
decreases the polarization. For example, scattered photons 
emerging from a = 0" were all scattered hy 6 = 180° and 
are therefore unpolarized (11(180°) = 0). 

Our radiative transfer code does not account for gas 
clumping in the IGM (in our simulation the gas density is 
a smooth function of radius r). Gas clumping affects the 
predicted polarization: Lya photons that resonate with gas 
inside denser clumps scatter more than calculated in our 
model, and are therefore polarized at lower levels (see FigU)). 
On the other hand, photons that resonate with gas in the in- 
terclump medium scatter less and are expected to me more 
strongly polarized. Since this interclump medium has the 
largest volume filling factor, gas clumping results in net 
higher level of polarization. The impact of gas clumping on 
the predicted polarization properties - at a fixed level of the 
ionizing background - can therefore be mimicked by increas- 
ing the level of the ionizing background somewhat. 

In summary, the polarization of resonantly scattered 
Lya photons in the IGM is low compared to the values de- 
rived in § 13.11 with the maximum polarization level reaching 
values of Pmax ;$7%. The polarization amplitude depends 
strongly on the local level of the ionizing background. This is 
because the polarization decreases with increasing effective 
optical depth, Tcb, which in turn increases with a decreasing 
local photoionization rate. 



3.3 Scattering in Optically Thick Clouds 

3. 3. 1 Description of the Model 

iDiikstra et al.l (|2006l ) performed Monte-Carlo calculations of 



0.6 



A 



0.4 



0.2 



V 




-0.2 



Figure 8. The level of anisotropy of the Lya radiation field, quan- 
tified by (n ■ er) (see text), is shown. For an isotropic radiation 
field (n ■ e^) = 0. Clearly, (n ■ e^) > and increases with radius. 
Photons that are scattered towards the observer were therefore 
preferentially moving radially outward prior to their last scatter- 
ing event. This Figure shows clearly that the Lya radiation field 
becomes increasingly anisotropic with radius, which results in an 
increasingly large fractional polarization with impact parameter 

(Figig. 



Lya radiative transfer through optically thick, spherically 
symmetric, collapsing gas clouds, which represent simplified 
models of protogalaxies in the process of their assembly. 
Here we compute the Lya polarizat ion properties for on e 
of the fiducial models presented by IDiikstra et al.l ()2006l i , 
in which a dark matter halo of mass Mtot = 5.2 x 10^^ M© 
collapses at 2 = 3 while it continuously emits Lya over a 
spatially extended region (see 'model 1.' of Dijkstra et al 
2006). 



3.3.2 Results 

The left panel of Figure [7] shows that the Lya emerges with 
a systematic blueshift, which is due to energy transfer from 
the gas to the photons, combined with a reduced escape 
pro bability for photons in the red wing of the line profile 
(see lDiikstra et al]|2006l . for a detailed discussion). The right 
panels of Figure [7] show the surface brightness profile (top), 
and the fractional polarization P(a) [bottom) as a function 
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of a. The fractional polarization again reaches values as high 
as Pmax = 35%. This may be surprising given the fact that 
each photon scatters numerous times before emerging from 
the cloud. As mentioned in § 12.11 Lya photons generally es- 
cape from an optically thick medium in the line wing, where 
they acquire larger polarizations than in the line core. Fur- 
thermore, each spherical shell inside the cloud must see a 
net outward flow of Lya photons, or else the cloud would 
not cool. 

This is quantified in Figure [8] which shows the radial 
dependence of the average value of n • e^. Here n denotes a 
photon's propagation direction prior to the last scattering 
event (as in Fig and er denotes a unit- vector that is 
pointing radially outward from the center of the cloud. In 
this convention n • er = 1 (n • er = — 1) for a photon that is 
moving radially outward (inward). For an isotropic radiation 
field (n • er) = 0. Clearly, (n ■ er) is positive at all radii, 
and increases with radius. In other words, photons that are 
scattered towards the observer were preferentially moving 
radially outward prior to the scattering event. Furthermore, 
this anisotropy in the radiation fielcQ increases with radius, 
which results in an increasingly large fractional polarization 
with impact parameter (Fig[7|. 

In Figure |5] we show the frequency dependence of the 
polarization at 3 different angular separations: a 0, 4.0, 
and 8.0 arcseconds. The thick solid line is the spectrum 
shown previously in Figure [T] indicating that the bulk of 
the photons emerge with frequencies — 10 < a; < 80. Within 
this frequency range, the polarization increases towards the 
blue. The reason for this wavelength dependent polarization 
is that Lya photons that are far in the blue wing of the line 
profile, appear even farther in the wing in the frame of the 
infalling gas. Therefore, the escape probability of Lya pho- 
tons increases, and the radiation field becomes increasingly 
anisotropic towards higher photon frequencies. This in turn 
increases the polarization towards higher frequencies. Lastly, 
it should be pointed out that the frequency dependence of 
the polarization applies mostly to the case a ~ 8.0 arcsec- 
onds, and similarly to the case a ~ ac that was shown in 
Figure [5] Hence, the frequency dependence of the polariza- 
tion is evident only where the overall (frequency averaged) 
polarization is significant (V ^10%). 



4 DISCUSSION 

4.1 Constraints on EoR and High-z Galaxy 
Formation 

The Lya emission needs to be spatially resolved in order 
for its polarization to be measured. Over the past decade, 
a large number of spatially extended Lya emitters, also 
known as Lya 'blobs', were discovere d jSteidel et al.|[200ol : 
iMatsuda etai]|2004l : [Saito et al.ll2"006l ). In fact, the size dis- 
tribution of Lya sources is continuous and Lya blobs simply 

■* We point out that the actual Lya radiation field that is 'seen' 
by each atom in the cloud is dominated by photons that are res- 
onantly scattering, and is almost completely isotropic. However, 
these photons have a negligible probability of escaping towards 
the observer and do not affect the observable Ly« polarization 
properties of the cloud. 
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Figure 9. The frequency dependence of the polarization at three 
different angular separations from the galaxy: o ~ 0, 4.0, and 8.0 
arcseconds. The polarization increases with separation (as was al- 
ready demonstrated in Fig[7ll. The thick solid line shows the spec- 
trum (arbitrary units, previously shown in Figure O. The bulk 
of the photons emerge with frequencies — 10 < a; < 80. Within 
this frequency range, the polarization increases towards higher 
frequencies, in sharp contrast with the outflow model discussed 
in S l3Tl 

repr esent the la rgest and rarest objects in this distribution 
( Mat suda et al . 2004). This implies that there is no shortage 
of spatially extended Lya sources in the Universe. 

As was shown in § [3l the polarization properties of spa- 
tially resolved Lya emission encode information about the 
medium through which the Lya photons propagate, espe- 
cially in combination with spectral information (Fig [5] and 
Fig [9]). For example, large fractional polarizations (tens of 
percent) are expected for Rayleigh scattering. An increase 
in the polarization towards longer wavelengths would be in- 
dicative of backscattering off an expanding outfiow. Such 
an inference could be particularly important because (i) it 
provides evidence that the extended Lya emission results 
from scattering (ii) it allows to better gauge the impact of 
resonant scattering in the IGM on the observed Lya line pro- 
file, since back-scattered Lya radiation is much less prone 
to resonant scattering in the IGM. Understanding the im- 
pact of resonant scattering on the Lya line profiles is essen- 
tial for (t ig) deriving constrain ts on the ionization state of 
the IGM (|Diikstra et al.ll2007al '). and ( iib) interpreting Lya 
emitters with large equivalent widths (|Malhotra fc Rhoadsl 
l2002l : lDiikstra fc Wvith3l2007l ). 

In this paper we have shown numerical examples for 
the polarization properties of Lya radiation that is either 
(i) resonantly scattered in the IGM; (ii) backscattered in 
an outfiow; or (Hi) associated with cooling radiation from a 
collapsing protogalaxy. There are other classes of spatially 
extended Lya sources that have not been discussed: 

• Gas in close proximity to a bright source of ionizing ra- 
diati on (e.g. a quasar) may 'glow ' in fluorescent Lya emis- 
sion i Haiman fc Rocs 2001; Wcidi nger et al.ll2004h . If the re- 
combining Lya emitting gas is sufficiently close to the quasar 
or the quasar is particularly bright, then this gas is optically 
thin to Lya photons and the emergent Lya radiation is not 
expected to be polarized. However, fluorescent Lya emission 
is also expected from Lyman limit systems and Damped Lya 
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systems (iGould fc Weinberd Il996l: ICantalupo et all l2005l : 
lAdelberger et alj l2006l : ICantalupo et alj l200'it l~ which are 
optically thick to Lya. Especially in cases where the gas 
is illuminated anisotropically, scattering could lead to a de- 
tectable polarization. 

• Spatially extended Lya emission may be generated 
by gas that was shock heated by a supcrwind-driven out- 
flow (e.g. .Taniguchi & Shiova 2000; Mori ot al. 200j). The 
backscattering mechanism described in ij |3.1l mav also oper- 
ate here, and this could lead to a measurable polarization 
amplitude. 

• Extended luminous Lya emission is common around ra- 
dio galaxies, a nd is thought to be po wered by a jet-IGM in- 
teraction (e.g. [chambers et al.lll99Cll ). Little is known about 
the source of the Lya emission and to what degree the Lya is 
scattered by residual neutral gas. Polarimetry may provide 
new insights on the nature of these sources. 

We have focused on spatially extended Lya sources, as 
the polarization averages out to zero for a symmetric unre - 
solved source. However, as was shown bv lLee fc Ahiil (|l998l ). 
anisotropic outflows (e.g biconical outflows) may still result 
in observable polarization, even if the Lya source is com- 
pletely unresolved. Hence, a measurable polarization may 
not be restricted to spatially resolved Lya sources. 

Lastly, we point out that accurate polarization mea- 
surements of faint high-redshift sources is difficult. It would 
therefore be prudent to obtain polarimetry of the brightest 
sources first. Polarization measurements of one of the bright- 
est known Lya Blobs (at z = 2.7, Dey et al. 2005) have al- 
ready been attempted (Prescott et al. 2008). This approach 
will likely lead to insights about the processes that regulate 
the transfer of Lya in and around high-redshift galaxies. 



4.2 Comparison to Polarization of Foregrounds 

The dominant foregrounds in the infrared sky (at wave- 
lengths A < 2^m), are R ayleigh scattered moonlight and 
starlight (e.g. iGlasd [T999I . p. 33), and are therefore polar- 
ized. However, line emitters can easily be extracted from this 
scattered foreground for which the spectrum is continuous. 
Here, we investigate whether Lya emitters can be distin- 
guished from foreground line emitters, such as [Oil] , [OHI] , 
Ha, and H/3 emitters, based on their polarization properties. 

In practice it has been possible to effectively discrimi- 
nate between these low-redshift interlopers and actual high- 
redshift Lya emitters base d on broad-band colors and on the 
shape of the Lya line (e.g. iKas hikawa et al ."2006''). However, 
the continuum spectrum of some Lya emitters with a large 
equivalent width is difficult to detect, and in some cases the 
Lya line is not asymmetric (see e.g the lower left panel of 
FigEJ. 

The possibility that high-redshift Lya emitters are 
highly polarized may provide another diagnostic that dis- 
tinguishes them from low-redshift interlopers. A low level 
of polarization {V <,2%) is expected for any extragalac- 
tic lin e emitter due to scatt ering by dust in the Milky- 
Way (jSchmidt fc Milled lloM ). Furthermore, spatially re- 
solved [OHI] line emission at higher levels of polarization 
{V ~ 4%) has been observ ed in a small fraction of Seyfert 
2 galaxies (|Goodrich|[l993 ). This polarization arises when 
the narrow line region is obscured. In this case, the only 



[OHI] photons that are observed are scattered towards the 
observer, either by dust or by electrons. This same mecha- 
nism could produce polarized [OH], Ha and H/3 line emis- 
sion. This suggests it is risky to distinguish high-redshift 
Lya emitters from low-redshift interlopers on the basis of 
polarimetry. On the other hand, the observed levels of po- 
larization in low-redshift interlopers are small compared to 
the values we find in our models, which suggests that high 
levels of polarization {V ~ 10%) are indicative of scattered 
Lya. 

Raman scattering of Ly/3 wing photons may produce 
highly polarized Ha emission in a radiative cascade of the 
form 3p ^ 2s ^ Is fe.g. lLee fc Yunlll998l : lYoo et al.ll2002l ). 
This mechanism requires the presence of large HI column 
densities A'hi ^lO'^" cm~'^, and has be en observed to oper- 
ate around symbiotic stars (jlkeda et al . 2004), a nd po ten- 
tially in close proximity to quasars I Lee fc Yunlll998l ). In 
distant galaxies, however, these emission regions would be 
completely unresolved, and the polarization averages out to 
zero. On the other hand, columns in excess of A^hi JjlO^" 
cm~'^ can easily exist in the neighborhood of galaxies, and it 
is not possible to rule out altogether the possibility of highly 
polarized Ha emission as a contaminant of polarized Lya 
emitters. Nevertheless, polarized Ha lines produced by this 
mechanism are typically broad (FWHM - 20 - 40 A) and 
symmetric, with wings exte nding into both the longer and 
shorter wavelengths (see e.g lYoo et al.|[200^ . their Fig 4-8); 
these features clearly distinguish them from Lya line emit- 
ters. 

Ground-based searches for high-redshift Lya emitters 
are complicated by the presence of bright OH emission 
lines, which are produced by the excited radical OH* in the 
Ionosphere, following the reaction H-I-O3 — > 0H*-|- 02 (e.g. 
Glass , 1999,') . These lines make up the Meinel bands l|Meinell 
1950al lbl)" and are typically much brighter than Lya sources. 
Furthermore, these lines exhibit large temporal and spatial 
variations in intensity d ue to the passage of density waves 
through the ionosphere jRamsav et al.lll992l ). However, the 
wavelengths of these lines are well known, and future instru- 
ments (such as e.g. the Dark Age Z Lya Explorer, Horton et 
al, 2004) intend to search for redshifted {z = 6.5 — 12) Lya 
lines between these lines using high-resolution (i? — 1000) 
spectrographs. 

Polarimetry may provide an interesting alternative 
method of suppressing OH-emission lin es. Although the OH- 
lines have been studied in detail (e.g iRamsav et al. I I1992I : 
Maihara ot al. 199^), we are unaware of any polarization 
measurements. If OH-lines are unpolarized, then polarime- 
try may be used to remove them efficiently 0. For example. 



^ The Meinel bands correspond to rotational and vibrational 
states of the OH-molecule. It is conceivable that the earth's mag- 
netic field introduces a preferred axis of rotation and/or vibration, 
which would induce some level of polarization to the lines. How- 
ever, any large-scale polarization of the foreground OH lines due 
to the Earth's magnetic field would not have structure on the 
scales of tens of arcseconds and would therefore be distinguish- 
able from the arc-like features of scattered Lya radiation around 
high-redshift galaxies. After removing any uniform polarization 
component, the remaining foreground should be unpolarized on 
the small angular scales of interest here. Identifying polarized Lya 
lines could in principle be optimized by using special purpose fil- 
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consider an unpolarized OH sky line that is detected at the 
N-a level after an integration time tint- Suppose we repeat 
the same observation, but we create a first image by passing 
the incoming radiation through a linear polarizer for 0. Stint • 
Next, we create a second image by observing the remaining 
0. Stint through a linear polarizer that is rotated by 90 de- 
grees compared to the first observation. If one subtracts the 
two images, then any unpolarized emission is removed down 
the level of ~ cr. 



5 CONCLUSIONS 

The high-redshift Universe is known to contain luminous 
Lya emitting sources such as gala^xies and quasars. The Lya 
photons that are emitted by these sources are typically scat- 
tered both in the interstellar medium and in the intergalac- 
tic medium. In this paper we have calculated polarization 
properties of this scattered Lya radiation. 

We used a Monte-Carlo Lya radiative transfer code, 
and endowed each Lya photon with a polarization vec- 
tor which is perpendicular to the photon's wave vector. In 
this formulation the Stokes parameters resulted from bin- 
ning together multiple independent photons. We differenti- 
ated between resonant and wing scattering, as these are de- 
scribed by different scattering matrices (see Appendix IX)) . 
Wing scattering is described by classical Rayleigh scatter- 
ing and in this case we advanced the photon's polarization 
vector in a scattering event following the well-tested ap- 
proach of Rybicki & Loeb (1999). However, resonant scat- 
tering is described by a superposition of Rayleigh scattering 
and isotropic scattering with corresponding weights of 1/3 
and 2/3 (§ 12.211. In this cas e , min or modifications to the 
approach of Rybicki fc Loebl ( 19991 ) are required when ad- 
vancing the photon's polarization vector (ij I2.3|) . 

We have applied our code to three classes of models 
which represent the diverse sets of environments around 
high-redshift galaxies: 



• First, we computed the polarization of the back- 
scattered LyQ radiation observed from galaxies surrounded 
by a superwind-driven outfiow. We have found that the frac- 
tional polarization may reach values as high as "Pmax ~ 40%, 
where the maximum polarization depends on parameters 
such as the speed of the outflow, ?;exp, and the column den- 
sity of neutral hydrogen atoms, A^hi (§ I3.1|l . In this case 
we have found the polarization to increase towards longer 
wavelengths. 

• Second, we considered resonant scattering in the inter- 
galactic medium (IGM) after reionization. Residual inter- 
galactic hydrogen can scatter up to ^90% of the Lya emit- 
ted by galaxies out of the line of sight. (This high fraction 
does not apply to galaxies with superwind-driven outflows 
where the Lya photons are scattered out of resonance before 
they encounter the IGM.) We have found the polarization 
of the scattered Lya halo around galaxies to be polarized at 
lower levels (Pmax ^7%) than the halos around galaxies with 
outflows. This follows from the fact that Lya photons scatter 
multiple times when they enter resonance, which isotropizes 



ters that search for the circular polarization pattern shown in 
Figure [2] around the sources of interest. 



the local Lya radiation field and reduces its net polariza- 
tion. The polarization decreases with the mean number of 
scattering events the photons experience (i.e. the polariza- 
tion decreases with an increasing effective optical depth in 
the Lya line) . For this reason, we expect the polarization of 
resonantly scattered Lya photons in the IGM to be largest 
around low redshift sources or around sources where the 
IGM is more highly ionized than average (e.g. in the vicin- 
ity of a bright quasar, see § 13. 2p . 

The overall lower levels of polarization of resonantly scat- 
tered Lya radiation in a reionized IGM are significantly 
lower than that of pre-reionization L ya halos aro und galax- 
ies embedded in a fully neutral IGM fcoeb fc Rvbicki 1999i ). 
This suggests that the polarization properties of Lya halos 
may help constrain the neutral fraction of the IGM during 
the epoch of reionization. 

• Third, we considered neutral collapsing protogalaxies 
that are emitting Lya cooling radiation. Despite the fact 
that these clouds are extremely optically thick to Lya pho- 
tons, and photons typically scatter ~ r times, we found the 
polarization to linearly increase toward the edge of the cloud, 
reaching a maximum amplitude of Pmax ~ 35%. (§[2I3J- The 
resulting polarization increases towards shorter wavelengths 
in contrast to the trend found in outflow models. 

Our results indicate that resolved high-redshift Lya emis- 
sion may be highly polarized under a variety of likely cir- 
cumstances. High-redshift Lya emitters are usually distin- 
guished from low-redshift line emitters on the basis of their 
broad-band colors and their asymmetric spectral line shape. 
The work presented here implies that polarization may pro- 
vide an additional diagnostic. Moreover, polarimetry has the 
potential to better remove the glow of infrared lines in the 
Earth's atmosphere, which would improve the sensitivity 
of ground-based observations to high-redshift Lya emitting 
galaxies outside the currently available redshift windows. 

The polarization properties of Lya radiation encode in- 
formation about the distribution and kinematics of neutral 
gas in and around galaxies. Polarimetry therefore comple- 
ments the constraints that are derivable from spectroscopy. 
Hence, from a theorist's perspective it is well worth to in- 
clude polarization in Monte-Carlo calculations of Lya radia- 
tive transfer. 

The redshifted Lya line has provided us with an im- 
portant window into the high-redshift Universe. This work 
suggests that in order to fully exploit the observations, one 
should focus on both Stokes parameters I and Q, rather 
than just /. Polarization measurements to better than 1% 
accuracy can be carried out by existing facilities such as the 
FOcal Reducer/low dispersion Spectrograph (FORSljfl on 
the VLT, the LRIS imaging spectropolarime ter at the W.M. 
Keck Observatory llCoodrich fc Cohenl|2 003'). and the CI AO 
polarimeter on Subaru ( Tamura et al.ll2003il . The results of 
this paper imply that Lya emitting sources would provide 
excellent future targets for these instruments. 
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APPENDIX A: POLARIZATION OF 
RESONANTLY SCATTERED LYq 

Al Lya Scattering Matrix 

The scattering matrix for a scattering process that is a su- 
perposition of Rayleigh and isotropic scattering of weight Ei 
and E2 may be written as (Chandrasckhar 1960, Eq 250-258) 



3 / cos^ 9 
^ - 2^M 1 



where 



7? 



1 1 
1 1 



(Al) 



(A2) 



Here, cos 6 = n ■ n' (§ 12. 3|) . and // and are the com- 
ponents of the scattered intensity parallel and perpendicu- 
lar to the plane of scattering, respectively. If the incoming 
light is assumed to be unpolarized, then Ii — Ir — and 
the phase function and degree of polarization are given by 
IChandrasckhar 1960, Eqs. 250-258) 



p{e) = E2 + ^Ei + ^Ei cos^ e, n{e) = 



i + lt 



(A3) 

Comparison with Eq. (|3} reveals that Lya resonant scatter- 
ing of unpolarized light corresponds to the case Ei — ^ and 
E2 = |. 

As was mentioned in § 12.21 the non-zero degree of po- 
larization arises fuUjQ from scattering between the levels 
l'S'1/2 — 2P3/2 — > lS'1/2. This transition is denoted by H 
(while K denote s the 15*1/2 — * '2P1/2 — > lS'1/2 transition). 
iHamiltonl l|l947f ) showed that for the H and K transitions 
{Ei,E2)h = (|, |) and (Ei,E2)k = (0, 1). Hence, scattering 



A simple way to see why scattering through the 2P1/2 level 
does not result in any polarization is that the angular part of 
the wavefunc tion describi ng the electron in the 2Pi/2-state is a 
constant (e.g. IWhitelll934 p. 134). Hence an atom that is excited 
into the 2Pi/2-state has 'lost' any memory of the direction and/or 
polarization of the photon that excited the atom. The photon that 
is emitted when the atom decays back to the ground state has no 
preferred direction (i.e. it is emitted with an equal probability in 
all directions). 




1215.65 1215.66 1215.67 1215 68 1215.69 
Rest frame X (A) 

Figure Al. The impact of quantum interference on the polar- 
ization in resonant Lya scattering. The frequency depe ndence of 
the parameter Ei (see Eo. lAll l. as calculated by Istenflgl (jl980|), is 
shown. One finds Ei = i and Ei = at the H and K resonance 
frequencies as expected. At sufliciently large separations from 
both resonances, Pi 1, which corresponds to pure Rayleigh 
scattering. Interestingly, Pi < for a range of frequencies be- 
tween the H and K resonance frequency. In practice, however, 
scattering at these frequencies does not occur often enough to 
leave an observable imprint (see text, and Fig IA2| |. The upper 
horizontal axis shows frequency in normalized units x (we arbi- 
trarily set X = a.t u = uk). 



by 90° via the H-transition results in n(90°) = |, as was al- 
ready mentioned in ii l2.2l Because scattering is twice as likely 
to occur via the H transition, the scattering averaged value 
becomes {Ei,E2) = i(Pi,£;2)K + |(£;i,P2)h = (§, |). 

Lastly, we point out that Rayleigh and isotropic scatter- 
ing by 90° results in 100 and 0% polarization, respectively. 
Since resonant scattering may be viewed as a superposition 
of Rayleigh scattering and isotropic scattering with corre- 
sponding weights of 1/3 and 2/3, one might expect 90° res- 
onant scattering to result in a fractional polarization of i, 
while the actual number is n(90°) = ^ < |- The main 
reason for this difference is that isotropic and Rayleigh scat- 
tering have different phase functions, and that i/ scattering 
by 90° occurs, then the probability that this was due to 
isotropic scattering is larger than |. Quantitatively, the an- 
gular dependence of the fractional contribution of Rayleigh 
scattering to the phase- function is given by Pr{9) = |(f + 
6')/(t| + 3 cos^ 6') = (1 + cos^ S)/{f + cos^ which 



4 V 12 

was used in i 



[23] 



A2 Quantum Interference as in Stenflo (1980) 



In § 12.21 the phase function and the degree of polarization 
caused by resonant and wing scattering were give n. Here 
we discuss the transition between these two regimes. IStenflol 
(1980) has calculated the value of Pi as a function of fre- 
quency (but note that Pi = W2 in the notation of Sten- 
flo 1980, Eq. 3.25). We have plotted Pi as a function of 
wavelength (bottom label) and frequency (top label) in Fig- 
ure [All 

Figure lAll shows that Pi = i at the H resonance fre- 
quency and Pi = at the K resonance frequency, which 
correspond exactly to the values quoted above. The Figure 
also shows that at sufficiently large separations from both 
resonances. Pi — > 1. In dimensionless frequency units, this 
asymptotic value is practically reached when \xat \ ^0.2. This 
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was the motivation for using the threshold x^it = 0.2 to sep- 
arate core from wing scattering in § 12.31 

Interestingly, i5i < for a range of frequencies between 
the H and K resonance frequencies. However, scattering at 
these frequencies does not occur often enough to leave an 
observable imprint. The reason for this is simple: the natural 
width of the line for both the H and K transitions is much 
smaller than their separation, i.e 7h,k ~ 10* Hz <^ uh — 
vk = 1.1 X 10^" Hz. Since the absorption cross-section in the 
atom's rest-frame scales as, a{v) oc [(!^ — t'H,K)^+7H,K/4]~^, 
a Lya photon is much more likely absorbed by an atom for 
which the photon appears exactly at resonance, than by an 
atom for which the photon has an energy corresponding to 
a negative W2. Quantitatively, the Maxwellian probability 
P that a photon of frequency x is scattered at frequency lat 
in the atom's rest-frame is given by 



P{x!,t,x)dx^t = 



-(aiat— a;) 



■KH{a, x) x'^f. + 



-dx. 



(A4) 



where H{a,x) is the Voigt function {H{a,x) = Tx/tq, see 
Eq[2l). This probability is shown in Figure IX2I foiFI a; = 3.3 
{solid line) and x = —5.0 {dashed line). Figure [X2I shows 
that for X — 3.3, photons are scattered either when they 
are exactly at resonance or when they appear ~ 3 Doppler 
widths away from resonance. The inset of Figure IX2I zooms 
on the region near a;at ~ 0.0. Clearly, if a photon is reso- 
nantly scattered, then it is more likely to be scattered when 
it is exactly at resonance then when it is, say, 0.05 Doppler 
widths away from resonance. For frequencies > 3.3 there 
are not enough atoms moving at velocities such that the 
Lya photon appears at resonance in the frame of the atom. 
Instead, the photon is scattered while it is in the wing of 
the absorption profile. This is illustrated by the dashed line 
which shows the case x — —5.0, for which resonant scatter- 
ing is less likely by orders of magnitude. 

From the above we conclude that (i) when a photon 
is scattered in the wing of the line profile, then Ei — 1 
and {a) when a photon is scattered resonantly, then it is 
scattered almost exactly at resonance, and therefore i5i = | 
(the scattering-averaged value). 



A3 A Note on the Polarization Dependence of 
the Scattering Phase-Function 

In this section we show that the procedure for propagating 
the wave and polarizations vector in a single Rayleigh scat- 
tering event that was outlined in § 12.31 naturally accounts 
for the polarization dependence of p{d) and n(&). 

We introduce a rectangular coordinate system in which 
a Lya photon propagates along the i— axis prior to scatter- 
ing, and in the x — z plane after scattering (see Fig IA3|I . In 
this coordinate system, n = (1,0,0), n' = (cos 61, 0, sin 6*), 
and e = (0, sin/3, cos/3), where the notation of § 12.31 was 
used. Here, /3 denotes the angle between the polarization 
vector e and the x — z plane, and 'i' denotes the angle be- 
tween e and n'. The components of the intensity parallel and 
perpendicular to the plane of scattering -prior to scattering- 
are given by /i — /cos^/3 and Jr — /sin^/3, respectively. 



° Note that the transition from core to wing scattering occurs at 
X ~ 3.3, see §0 
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Figure A2. The probability that a photon of frequency x is 
scattered by an atom such that it appears at a frequency x^t 
in the frame of the atom. The solid and dashed lines correspond 
to a; = 3.3 and x = —5.0 respectively, P(xat,2^) (Eg. IA4T l. For 
X = 3.3, photons are either scattered by atoms to which they 
appear exactly at resonance (see inset), or to which they appear 
~ 3 Doppler widths away. For x = —5, resonant scattering is less 
important by orders of magnitude. In combination with Fig lAll 
this figure shows that if a photon is resonantly scattered then Ei 
is either or |. 




Figure A3. The coordinate system that is used to investigate the 
dependence of the scattering phase function on the polarization 
of the incoming radiation. 



From the scattering matrix it follows that the phase func- 
tion for this polarized radiation is given by 

p{e)dn = ^ cos^ e cos^ I3 + ^ sln^ f3j dQ. (A5) 

Using that cos ^ = n' ■ e = sin 9 cos /3 we obtain 



p(*)df7 = ^ sin^ 



(A6) 



Hence, the probability that a photon is emitted into solid 
angle dQ,{6, (j>) scales simply as sin^ \E' (where ^ — 'if{6, 4>)). 
This phase-function was used in this paper for Rayleigh scat- 
tering (sin^ ^ = 1 — {e ■ n')^). Hence, this method naturally 
incorporates the polarization dependence of scattering phase 
function. Furthermore, the intensity parallel and perpendic- 
ular to the plane of scattering after scattering are given by 
// = I /cos /3cos 9 and = |/sin /3, respectively. The 
polarization vector of this scattered radiation is therefore 
indeed given by e', which was obtained by projecting e onto 
the plane perpendicular to n' (see §[2]3}. 



iT^ 9nn7 RA« MNTRAS nnn miTZl 



